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Abstract 

During the 2008 run RHIC provided high luminosity in both p+p and d+Au 
collisions at ^s^n = 200 GeV. Electromagnetic calorimeter acceptance in 
STAR was enhanced by the new Forward Meson Spectrometer (FMS), and 
is now almost contiguous from —1 < rj < 4 over the full azimuth. This large 
acceptance provides sensitivity to the gluon density in the nucleus down 
to x ~ 10~ 3 , as expected for 2 — y 2 parton scattering. Measurements of the 
azimuthal correlation between a forward 7r° and an associated particle at large 
rapidity are sensitive to the low-x gluon density. Data exhibit the qualitative 
features expected from gluon saturation. A comparison to calculations using 
the Color Glass Condensate (CGC) model is presented. 
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1. Introduction 

Two particle azimuthal correlations represent a powerful tool for charac- 
terizing the transitional region between dilute and saturated partonic sys- 
tems. Inclusive particle production measurements from early RHIC runs 
show general agreement with pertubative QCD for p+p collisions at a/s=200 
GeV over a broad range of rapidity. Results from d+Au collisions show 
Cronin type enhancements at midrapidity and strong suppression for for- 
ward particle production jl|, 0, 0], suggestive of parton saturation, although 
consistent with other explanations. The broad range of longitudinal momen- 
tum fraction x that is averaged over in inclusive production is more selec- 
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tively probed in two particle correlations, especially when the pseudorapidity 
(77 = —log(tan(9/2))) of one of the two particles is varied. 

In pQCD at leading order, particle production in high energy hadronic 
interactions results from the elastic scattering of two partons (2—7-2 scat- 
tering). Each initial state parton carries a fraction of the colliding hadron 
momentum, as given by universal parton distribution functions. The scat- 
tered partons fragment to the observed hadrons. Although complexities arise 
in calculating two-particle correlations using pQCD, 2—7-2 partonic scatter- 
ing leads to the expectation of events with back-to-back jets. When high-pr 
hadrons are used as jet surrogates, we expect the azimuthal correlations of 
hadron pairs to show a peak at A0 ~ it from back-to-back jets, and a peak 
at A0 ~ when At] between the two particles is smaller than a typical jet 
size. Transverse momentum effects will broaden these peaks in the azimuthal 
correlation function 

When the parton density increases, the basic dynamics for the particle 
production can change. Instead of elastic 2 — > 2 scattering, the particle pro- 
duction can proceed by the interaction of a probe parton from the deuteron 
beam with multiple gluons from the Au beam. At sufficiently high parton 
densities, the transverse momentum from the fragments of the probe parton 
may be compensated by several gluons with lower p?. Two particle azimuthal 
correlations are expected to show a broadening of the back-to-back peak (loss 
of correlation: 2 — > many processes) and eventually to disappear (monojet). 

Different models try to describe hadron production from dense (nuclear) 
targets by including non-linear contributions. One approach |4j extends per- 
turbative QCD factorization by adding contributions of quarks scattering 
coherently off multiple partons in the target. In the Color Glass Condensate 
(CGC) framework 0,0, 0,S, [9}, the hadronic wave-function is saturated as a 
consequence of gluon recombination. At very low values of the longitudinal 
momentum fraction x of the probed gluons the occupation numbers become 
large. In both approaches the probe scatters coherently off the dense gluon 
field of the target, which recoils collectively, leading to a modification in Atp. 

2. Experimental Setup 

The STAR detector at RHIC provides the suitable configuration for study- 
ing two particle azimuthal correlations. Electromagnetic calorimeters pro- 
vide nearly hermetic coverage for 2n in azimuth and —1.0 < rj < 4.0. The 
STAR calorimetry includes the Barrel ElectroMagnetic Calorimeter (BEMC: 
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Figure 1: On the left: schematic view of the STAR experimental hall. Gold beam is 
coming from the West (right in figure) side. On the right: PYTHIA [10] simulation of 
di-pion production at large 77 p+p collisions at y/s = 200 GeV. The r\ of the associated 
particle is strongly correlated to the x value of the soft parton probed in the partonic 
scattering. Figure from [111 ]. 

\rj\ < 1.0), the End-cap ElectroMagnetic Calorimeter (EEMC: 1.0 < 77 < 2.0) 
and the newly installed Forward Meson Spectrometer (FMS: 2.5 < rj < 4.0). 
These, together with the STAR Time Projection Chamber (TPC), allow 
measurements of correlations of different particles over a broad Arj x Atp 
range. 

Signatures of possible saturation effects can be investigated by compar- 
ing azimuthal correlations in p+p and d+Au interactions. The first case is 
used as a reference: the colliding systems are dilute and no saturation effects 
are expected at these energies. In d+Au, on the contrary, the dense nuclear 
matter is probed and azimuthal correlations are expected to broaden when 
the saturation scale is reached. No large final state rescattering effects (as 
in Au-Au collisions) are expected in d+Au. The measurement of di-hadron 
azimuthal correlation functions involving a forward neutral pion, and vari- 
able 7] for the second hadron, are expected to reveal signatures of saturation 
when the gluons are at sufficiently low x. The STAR trigger systems allow 
us to select events with a leading ir° within the FMS acceptance. In a pQCD 
description, neutral pions are produced in the forward direction by the inter- 
action between a large-a; parton (most likely a valence quark) from the beam 
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heading towards the FMS (either p or d) with a \ow-x gluon from the other 
beam (either p or Au). The correlation probability is then calculated by 
detecting a second particle in any of the STAR sub-detectors. Since the sat- 
uration scale depends both on the momentum fraction x of the probed gluon 
and on the scale Q 2 , one can approach the saturation region by requiring 
an associated particle with increasing pseudo-rapidity, or by selecting events 
characterized by a lower p? respectively. Additional insights can be achieved 
by selecting correlations based on the centrality of the collision, since a larger 
effect is expected when the denser part of the nucleus is probed. 



3. Results 



A plan of azimuthal correlation measurements has been systemati 



cally pursued at STAR. The objective has been to determine if the boundaries 



of the saturation region are accessible at RHIC energies [12| and establish 
the effects on particle production. The pt dependences of high density nu- 
clear effects have been studied using prototype forward calorimeters (Forward 
Pion Detector, FPD/FPD++, [1]) and recently the FMS Comparison 
of Aip between a forward 7r° and a mid-rapidity hadron in p+p and d+Au 
shows a significant broadening in the back-to-back peak in d+Au. Such ef- 
fect appears to be stronger as the pt of the particles decreases, as expected 
from saturation models. The x dependence is accessed by comparing the 
effects on A(f while varying the pseudo-rapidity of the associated particle. 
Correlations between two forward pions lij show indeed a much stronger 



broadening than in forward + mid-rapidity correlations. Measurements in 
the intermediate rj region (using the STAR EEMC) are ongoing. The cen- 
trality dependence of the effect was studied for all 77 and px regimes. In this 
analysis, the impact parameter b of the collision is evaluated in relation to 
the multiplicity, measured as the sum of charges (J2Qbbc) in the event by 
the east side of the STAR Beam-Beam Counter (BBC), facing the Au beam. 
The multiplicity classes for peripheral (0 < J^Qbbc < 500) and central 
d+Au collisions (2000 < J^ Qbbc < 4000) correspond to (b) = (6.8 ± 1.7) 
fm and (b) = (2.7 ± 1.3) fm respectively, as indicated from a HIJING 1.383 
fl5| simulation. 

Figure [2] shows the (efficiency uncorrected) probability to find an asso- 
ciated pion given a trigger 7r°, both in the forward (FMS) region. Neutral 
pions are reconstructed from pairs of photon clusters within the FMS fidu- 
cial volume. Pion candidates are required to have an invariant mass in the 
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Figure 2: Uncorrected coincidence probability versus azimuthal angle difference between 
two forward neutral pions in p+p collisions (left) compared to peripheral (center) and cen- 
tral d+Au collisions (right). Data are shown with statistical errors and fit with a constant 
plus two Gaussian functions (in red). CGC expectations fli| have been superimposed (in 
blue) to data for central d+Au collisions. 



interval 0.05 < M 77 < 0.25 GeV/c 2 . The pair with the largest pt is selected 
as the leading (trigger) pion and its azimuthal coordinate is compared inclu- 
sively with those of all other (associated) pion candidates. Trigger particles 
are selected with a transverse momentum of p^ 9 ^ > 2.0 GeV/c while asso- 
ciated particles are required to have p^ 9 ^ > p^ sscS> > 1.0 GeV/c. In Figure 
[2] we compare A(p for p+p interactions with peripheral and central d+Au 
collisions. All the distributions present two signal components, surmount- 
ing a constant background representing the underlying event contribution 
(larger in d+Au). The peak centered at Aip = (near-side peak) represents 
the contribution from pairs of pions belonging to the same jet. It is not 
expected to be affected by saturation effects, therefore it is a useful tool to 
check the effective amount of broadening in the away-side peak. This second 
peak, centered at Aip = n, represents the back-to-back contribution to the 
coincidence probability which is expected to disappear in going from p+p to 
d+Au if saturation sets in. Data are fit with a constant plus two Gaussians 
centered at A(p = and Aip = n. Data show that the width of the near- 
side peak remains nearly unchanged from p+p to d+Au, and particularly 
from peripheral to central d+Au collisions. The away-side peak, instead, 
presents strong differences between peripheral and central d+Au collisions. 
Peripheral d+Au collisions show back-to-back contribution like in p+p, even 
though apparently smaller in height relative to the near-side peak than in 
p+p, but broadened. Central d+Au collisions show a substantially reduced 
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Figure 3: Comparison between parton density functions and their effect on azimuthal 
correlations. Left: comparison of density distributions of gluons using Duke-Owen (blue), 
CTEQ3L (red) and CTEQ5L (black) PDF sets. Right: uncorrected forward di-pion cor- 
relations in p+p as simulated in PYTHIA 6.222, using the three different PDF sets. 



away side peak that is significantly broadened. 

The rightmost panel of Figure |2] shows also a comparison (in blue) with 
theoretical expectations using the CGC framework. Inclusive particle pro- 
duction in interactions between a dilute system (deuteron) and a saturated 
target (Gold) has been calculated using a fixed saturation scale Qs [12| and, 
consequently, used to compute di-hadron correlations [l6| . Calculations con- 
sider valence quarks in the deuteron scattering off low-x gluons in the nucleus 
with impact parameter b = 0. CGC calculations show qualitative consistency 
with data in their expectations of a strong suppression of the away-side peak 
in central d+Au collisions. 



4. Systematics 

Further studies have been performed on this analysis with the goal of 
evaluating the expectations of "conventional" calculations of azimuthal cor- 
relations, using PYTHIA 6.222 and HIJING 1.383. Forward di-hadron cor- 
relations are expected to be sensitive to gluon density, initial-state showers 
and intrinsic transverse momentum fcy of the partons. HIJING models d+Au 
collisions by using a version of PYTHIA, and a Glauber model for the nuclear 
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Figure 4: Uncorrected azimuthal correlations between two forward neutral pions as sim- 
ulated in p+p by PYTHIA 6.222 (left). On the central and rightmost panel: PYTHIA 
events have been embedded into minimum bias d+Au data to emulate the effect of the 
larger multiplicity in d+Au. Each probability is scaled by a common arbitrary factor. 



effects. The latest available HIJING versions present only a limited choice of 
parton distributions functions (PDF), all of which predate HERA discovery 
of the rapid growth of the gluon distribution at low-x. A systematic study 
of the dependence of forward di-pion azimuthal correlations on the PDF has 
been made using PYTHIA 6.222. Figure |3] shows the comparison between 



three different sets of parton distributions: the Duke-Owen set 17 , which 



predates HERA data [18|, [19|] and therefore does not include a rise in the 
gluon density at low-x, and the two sets CTEQ3L and CTEQ5L [201, the 
latter with a steeper rise of the gluon distribution at low-x. The right panel 
in Figure [3] shows the azimuthal correlations expected from these three PDF 
sets in p+p interactions. Data appear to be consistent in both peak areas 
and widths with a set of PDF that presents a rapid rise of the gluon density 
at low-a;. 

The effect on azimuthal correlations of the larger multiplicity in d+Au 
compared to p+p events has been studied. For this purpose, simulated di- 
pion PYTHIA events have been embedded into minimum bias d+Au data. 
Embedding has been preferred to pure HIJING simulations for reasons dis- 
cussed before. Figure @] shows the results of these studies. Azimuthal cor- 
relations between two forward pions are simulated using PYTHIA only and 
compared to PYTHIA events embedded into minimum bias peripheral and 
central d+Au interactions. The effect of the additional multiplicity in d+Au 
affects particularly the reconstruction of the associated tt°, causing an higher 
combinatorial background which reflects in the higher background level seen 
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Figure 5: Off-peak anaiysis. On the left panel invariant mass spectra for leading and 
sub-leading pion candidates are shown. The off-mass selection for sub-leading pions is 
indicated in blue. Following, uncorrected azimuthal correlations between two forward 
pion candidates (center) and between a forward pion candidate and a pair of clusters with 
invariant mass in the range 0.25 < M 11 < 0.45 GeV/c 2 (right). 



in d+Au in simulations as well as in data. The features of the near-side peak 
are reproduced in both p+p and d+Au in this embedding simulation. On 
the other side, the large suppression of the away-side peak noted in data, in 
particular in central d+Au collisions, is not reproduced in this simulation. 
We conclude that the away-side peak suppression is not the result of the 
multiplicity increase in central d+Au collisions relative to p+p collisions. 

As a further test, the impact of the combinatorial background on the 
shape of the azimuthal correlations has been studied. To do so, a pair of 
clusters with an invariant mass in the range 0.25 < M 77 < 0.45 GeV/c 2 , 
relatively far from the nominal value of the ir° mass, has been used as asso- 
ciated particle. Figure |5] shows the azimuthal correlation between a leading 
71° and such associated pairs. With this mass selection, Aip appears to be 
quantitatively consistent, in the widths of both peaks, with the "on-peak" 
correlation. Such small impact is expected for jet-like correlations where the 
combinatorial background is produced between pairs in the same jet. 

5. Conclusions 

Two particle correlations are a more sensitive probe of the transition re- 
gion between dilute and saturated systems than inclusive particle production. 
Thanks to a rich d+Au RHIC run in 2008, STAR is pursuing its objective 
of mapping the boundaries of the saturation region. Two particle azimuthal 
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correlations have been systematically studied as a function of p?, T] and multi- 
plicity When two hadrons are reconstructed in the forward region, where the 
lowest x in the gluon distribution is probed, comparisons with simulation in- 
dicates that p+p data are sensitive to the low- a; gluon density. Measurements 
over p+p and d+Au data show a strong suppression of the away-side peak in 
central d+Au collision compared to p+p, qualitatively consistent with CGC 
calculations. Systematic studies show that such suppression is not expected 
to arise from additional multiplicity in d+Au and it is not significantly af- 
fected by the presence of combinatorial background in the invariant mass 
di-photon spectrum. 
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